Charge-ordered ferromagnetic phase in manganites 
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A mechanism for charge-ordered ferromagnetic phase in manganites is proposed. The mechanism 
is based on the double exchange in the presence of diagonal disorder. It is modeled by a combination 
of the Ising double-exchange and the Falicov-Kimball model. Within the dynamical mean-field 
theory the charge and spin correlation function are explicitely calculated. It is shown that the system 
exhibits two successive phase transitions. The first one is the ferromagnetic phase transition, and 
the second one is a charge ordering. As a result a charge-ordered ferromagnetic phase is stabilized 
at low temperature. 
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There has been much recent interest in the proper- 
ties of doped manganese oxides Ri_a;Aa;Mn03 (R=rare 
earth, A=Ca,Sr)i^ These materials present a very rich 
phase diagram involving phases with spin, charge and 
orbital order. The physically relevant electrons in man- 
ganites are those from the Mn 3cZ levels, which are split 
by the cubic crystal field into triply degenerate t2g levels 
and higher-energy doubly degenerate levels. Electrons 
from the Cg levels are able to hop between Mn sites and 
form a conduction band. Electrons from the t2g levels are 
localized. The itinerant electrons and local spins are cor- 
related by the double-exchange (DE) mechanism?2i4 The 
main feature of the DE is a cooperative effect where the 
motion of an itinerant electron favors the ferromagnetic 
(FM) ordering of local spins and, vice versa, the presence 
of the FM order facilitates the motion of the itinerant 
electron. The DE model qualitatively describes some of 
the magnetic properties of manganiteS(2i^ and provides 
a well-established starting point toward comprehensive 
understanding of the phase diagram of manganites. 

Recently experiments have shown that beside the FM 
order a charge order can exist in the mang anitesM The 
charge order exists in regions with no net magnetization 
and, surprisingly, can also occur in FM regions J. Doping 
of ions creates Mn''+ holes in a Mn'^+ background. 
The presence of two valence states Mn'^+ and Mn^'^ 
may lead the compounds to a charge-ordered (CO) state 
for appropriate doping. However, the DE model alone 
cannot explain the CO state which coexists in the FM 
phase. In principle, the nearest-neighbor Coulomb repul- 
sion may stabilize a CO state. However, a large nearest- 
neighbor repulsion likely destabilizes the homogeneous 
FM state and may produce a checkerboard charge or- 
der in three directions.^ Another possible mechanism 
for the CO phase stabilization is the coupling of itin- 
erant electrons to the Jahn- Teller distortions. However, 
the electron Jahn- Teller phonon coupling can only sta- 
bilize a CO-FM state where the CO phase transition 
occurs before the FM transition.^ At half filling exper- 
iments have only observed a charge order below the FM 
transition temperaturei- -^ Therefore the Jahn- Teller cou- 
pling is unlikely responsible for the appearance of the 
CO-FM state at least at half filling. In this paper we 



present a possible alternative explanation for the CO- 
FM state in the manganites. The key idea is an inter- 
play of the DE and randomness of the A-site substitu- 
tion. The randomness is inevitably introduced by ex- 
periments. The importance of the randomness has been 
discussed both experimentally and theoreticallyi^ The 
randomness can substantially decrease the critical tem- 
perature of the FM transitionifiiifiiii Here we will incor- 
porate the randomness of A-site substitution into the 
DE model. For simplicity, we adopt the randomness by 
A-site substitution as a random local potential of the 
itinerant electrons, although the randomness may cause 
other effects, for instance, randomness of the hoping or 
exchange integralii^ It is well known that the diagonal 
disorder with binary distribution can be modeled by the 
Fahcov-KimbaU (FK) modelJ^ Although the FK model 
is simple, it contains a rich variety of phases. In par- 
ticular, it illustrates the disorder-order phase transition 
driven by electron interaction^^ii^ Incorporating the di- 
agonal disorder of the FK type into the DE model, one 
may expect that a disorder-order phase transition could 
present. When the phase transition occurs, a CO-FM 
phase may be stabilized at low temperature. In order to 
detect the phase transition we study the charge and spin 
response of system by using the dynamical mean-field 
theory (DMFT).i6 The DMFT has extensively been used 
for investigating strongly correlated electron systemsii^ 
Within the DMFT we explicitely calculate the charge 
and spin correlation function. We find that the system 
stabilizes a CO-FM state at low temperature. 

The system which we study is described by the follow- 
ing Hamiltonian 

<'ij>,(T 1(7 i 

E^'^Wi + U'^niaWi, (1) 

i ia 

where c|g.(Cj^) is the creation (annihilation) operator of 
an itinerant electron with spin a at lattice site «; t/Vd is 
the hoping parameter of the itinerant electrons. Here we 
have rescaled the hoping parameter with the dimension 
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d of the system. is the z component of local spin at 
lattice site i, and for simplicity, it takes two values ±1. 

= ('^n ~ nii)/2, = cl^c^^, Wi is a classical vari- 
able that assumes the value 1(0) if site i is occupied (not 
occupied) by A ion. U is the disorder strength and is 
mapped onto the difference in the local potential which 
splits energetically favor of Mn^+ and Mn^+ ions. The 
expectation value x = is the number of 

lattice sites), corresponds to the concentration of unfa- 
vorable Mn^+ sites. The chemical potential /i controls 
the carrier doping, while controls the fraction of the 
sites having the additional local potential. We shall use 
the condition n + x = 1, where n = X]icr('^io-)/-^ is the 
electron doping. This condition determines for each 
doping n. The third term of Hamiltonian |^ is the Hund 
coupling of itinerant and local electrons. For simplicity 
we only take into account the Ising part of the Hund 
coupling. This simplification does not allow any spin-flip 
processes, which can be important at low temperature 
where spin-wave excitations may govern the thermody- 
namics of the system. However, in the DE processes the 
spin of itinerant electron ferromagnetically aligns with 
the local spin, hence, the Ising part of the Hund interac- 
tion plays a dominant role. The DMFT calculations for 
the DE model with classical local spins show that the sim- 
plification of the Hund coupling does not change the self 
energy of the single-particle Green functioni^i Moreover, 
within the DMFT the numerical results for quantum lo- 
cal spins do not show a significant difference from the 
ones for classical local spinsi^^ Thus, one expects that 
within the DMFT the simplification of the Hund inter- 
action does not result in a serious backwardness. Jh is 
the strength of the Hund coupling, and in the following 
we will take the limit Jh — > oo. The first three terms of 
Hamiltonian iQJ constitute a simplified DE model. This 
simplified model captures the most essential ingredient 
of the DE processes. The last two terms of Hamiltonian 
(0 describe a binary randomness of the A-site substitu- 
tion. They together with the hoping term form the FK 
model. It is well known that within the FK model the 
U term induces a disorder-order phase transition.i^ii^ At 
low temperature a checkerboard ordering phase is stabi- 
lized. Hence, the model may display an interplay of 
the FM and CO phase. 

We solve model (HJ by the DMFT. The DMFT is based 
on the infinite-dimension limit. In the infinite-dimension 
limit the self energy is pure local and has no momentum 
dependence. The Green function of itinerant electrons 
satisfies the Dyson equation 



Go-(k, iLj„) = 
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e(k) + fi- Yja{i(^n) ' 
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where a;„ = 7rT(2n -t- 1), e(k) = ~2t^^^-^ cos(fcj), and 
So-(iWn) is the self energy. In the infinite-dimension limit 
the bare density of states of itinerant electrons becomes 
p{e) = exp(— e^/t^)/-y7ri and we take t as the unit of 
energy {t = 1). The self energy is determined by solving 
an effective single-site problem. The effective action of 



this problem is 

SoS = J2j dTdT'cl{T){-^S{T - t') + A^T - t'))c4t') 

(7 

+ / dTc\{T)[ — + Uw — Jh(tS^)c„{t) + EyjW, 

where A„(t) describes the effective medium. This effec- 
tive single-site problem can exactly be solved. Indeed, 
the dynamics of the localized spin and impurity w 
involved in the effective action are independent, hence, 
we could independently take the trace over and w in 
calculating the partition function. This is similar to the 
DMFT solving of the FK modelpi^ We obtain the local 
Green function in the limit Jh oo 
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where Z^{iujn) ~ ii^n + /i — Ao-(iwn), and 
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exp — jiE^ (pi — a) 

Za'{iu)n) - a'U 
Z„{iujn) - cdJ 



with a = 0, 1. In taking the limit Jh ^ oo in deriving 
Eq. Q we must first renormalize the chemical potential 
/i ^ ^ -I- Jh . The self energy is determined by the Dyson 
equation for the effective single-site problem 
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Within the DMFT, the local Green function must coin- 
cide with the single-site Green function of the original 
lattice, i.e., 



^^G^(k,iw„). 



(5) 



Eqs. ijlll-ISIl form the complete set of equations, which 
self consistently determine the self energy and Green 
function. 

We are interested in calculating the charge (c) and spin 
(s) correlation function 

x'^'Hi,j) = {{Sn,^ ± Sn.a){Sn,^ ± Sn,^^)) (6) 

in the homogeneous paramagnetic (FM) phase {Srii^ — 
nia — {nia))- In order to calculate the charge and spin 
response of system one has to introduce an external field 
into the Hamiltonian. The charge and spin correlation 
function can be obtained by differentiating the Green 
function respected to the external field, and then tak- 
ing the zero limit of the fieldii^ Following the standard 
techniques^ one can express the correlation functions 
in the terms of charge (c) and spin (s) susceptibility 
T^c(s)j'q^ jj^^j in momentum space 



X^(«)(q) 



(7) 
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The charge and spin susceptibihty can be obtained by 
differentiation of the Green function.™ We obtain 



0=0,1 



(8) 



where xo(q,«w„) ^ J2]iC^i'^ + 'ii'^^n)G(k,iLi!n). The ma- 
trix r'^(*''(q) satisfies the fohowing equation 

B=(^)(q)r(^)(q)=g^(^)(q), (9) 

where I3'^^^\q) and Q'^^^-'(q) have the foUowing elements 



E 



A 



l^'\lLU„)(^GHliU^){ 



1 - G'^{iu;n){^^)^ + ri{(\,iujr,)G{iujn) 
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with a, (3 — 0,1. In deriving Eqs. lP)- ((TT)) we 

have used the standard conversion^^ [xo{(l,iuJn)]~^ — 
[Gcr{iuJn)]^'^ + ??(q, i'^n)['-'o-(*'^n)]~"'", and introduced 

quantity v4a'^''(ia;„) — SWa] /SZi{iujn)iSWa^ /6Zi{iujn)- 
In the infinite dimension Umit all of the wave vector de- 
pendence of xo(q: ii-^n) and 77(q, icUn) include in the term 
X(q) = Hence, the spin and charge cor- 

relation function only depend on momentum via ^(q)- 
Each of the derivatives appearing in Eqs. can di- 

rectly be calculated from the DMFT solution of Eqs. Q- 
0. In such the way 13'^^^'>{q) and Q'^^^^q) are calcula- 
ble once the self-consistent equations of the DMFT are 
solved. Equation Q reveals that r'^^''''(q) will diverge 
at a temperature where the determinant of ]3^^^^q) van- 
ishes, while Q'^^'^\ci) remains finite. This results in an 
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FIG. 1: The critical temperature Tc as a function of doping 
n for U = 0.5 (squares), U — 1 (circles). The filled (open) 
symbols are Tc of the FM (CO) phase transition. The solid 
line is Tc of the FM transition without disorder (U — 0). 
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FIG. 2: The critical temperature Tc as a function of U for 
various doping n. The filled (open) symbols are Tc of the FM 
(CO) phase transition. 



unphysical change of the sign of corresponding correla- 
tion function x'^^^'Hq) so that the assumption of the ho- 
mogeneous PM phase fails for lower temperature. By a 
similar way one could also calculate the spin correlation 
function xs(q) of local spins. After some calculations we 
obtain 



xs(q) = 7o(q) +7i(q), 

T 



^ S^^(q)7Mq: 

13=0,1 



(12) 
(13) 



From Eqs. lO, p2|) - (|13|) one can see that the spin corre- 
lation function of itinerant electrons and local spins will 
diverge at the same temperature where the determinant 
of i?^(q) vanishes. This means that the spin of itiner- 
ant electrons parallel aligns with local spin, and thus is 
an important feature of the DE. We calculate the charge 
and spin correlation function 10 by solving the DMFT 
set of self-consistent Eqs. ©-(jSl- We are only interested 
in the FM and checkerboard CO phase stability. Hence, 
we only calculate the spin correlation function at Xc^ — 1 
and the charge correlation function at = — 1. It is 
found that the spin correlation function x^i^i^ — 1) al- 
ways diverges at a critical temperature. This is the sig- 
nal of the FM phase transition. The charge correlation 
function x'^(^q — ^1) only diverges for [/ ^0. This 
means that without the disorder the system always is 
homogeneous. In Figs. ^ and |2 we present the critical 
temperature Tc of the FM and CO phase transition as a 
function of doping and disorder strength. At half filling 
n = 0.5 both critical temperatures reach their maximal 
value. The Tc of the FM transition always decreases with 
increasing disorder strength. This means that the disor- 
der substantially decrease Tc of the FM transition. ^■-'^"■li 
At the same time, with increasing U , Tc of the CO phase 
transition first increases, reaches its maximal value, and 
then decreases. The behavior of Tc of the CO phase tran- 
sition is similar to the one in the FK model. ii> At very 
strong disorder (U ^ 1) the two critical temperatures 
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FIG. 3: Upper panel: the temperature dependence of the 
spin magnetization tua {jtib) of sublattice A(B) (the dot and 
dashed line), and of the charge-order parameter A (the solid 
line) . 

Lower panel: the temperature dependence of the free energy 
F. The solid, long dashed, short dashed, and dot-dashed line 
are the free energy in the CO-FM, homogeneous FM, CO-PM, 
and homogeneous FM phase, respectively {U = 1, n = U /2). 

approach to a same value . One also notices that Tc of 
the CO phase transition always is smaller than the FM 
transition temperature. Thus, one may expect that the 
CO state is stabilized in the FM phase at low temper- 
ature. However, this CO phase stability is respected to 
the homogeneous PM phase, and for safety we also study 
an inhomogeneous phase. We divide the lattice into two 
penetrating sublattices A and B. This lattice division 
allows us to study the checkerboard CO phase. By using 
the standard techniquei^ the matrix Green function can 
be written in the following form 

G-i(k,zc.„) = ('*'"" + '^~,,^,^(*"") , ^') 
where Tia^'^^iuJn) is the self energy of the Green func- 



tion of sublattice A{B). The self energies are deter- 
mined by solving the effective problem of single site 
of the sublattices We find that at low temperature 
a checkerboard CO-FM state is stabilized. We plot 
the magnetization m^(s) — '^J2i£A{B)i^i) /-^ '^^ 
lattice A (B) as a function of temperature in Fig. 13 
(upper panel). In this figure we also plot the temper- 
ature dependence of the charge-order parameter A — 

(E»ga,<t("-*'t> - JljEB^ai'^i'^))/^- It shows that below 
a critical temperature the magnetizations of both sub- 
lattices exist. They equal to each other until another 
critical temperature, where the charge-order parameter 
exists. At low temperature the system is in the checker- 
board CO-FM state. In this phase the charge order coex- 
ists in the FM state, as experimentally observedpi In the 
way the system exhibits two successive phase transitions. 
Initially the system goes to the homogeneous FM phase, 
and after that to the checkerboard CO phase, that the 
CO-FM phase is stabilized. We also calculate the free en- 
ergy of the system. The free energy can only be expressed 
in terms of local quantitiesiifi We plot the temperature 
dependence of the free energy F in Fig 13 (lower panel). 
It shows that the CO-FM state has lowest free energy, 
hence the state must be stabilized at low temperature. 

In conclusions, we have proposed a mechanism for the 
CO-FM phase which has recently been observed. The 
mechanism is based on a combination of diagonal disor- 
der and a simple DE model with local Ising spins. Em- 
ploying the DMFT we have calculated the charge and 
spin correlation function. It is found that the FM and 
CO state are stabilized at low temperature. As a result 
the checkerboard charge order can occur in the FM state. 
However, the manganites are too complicated a system 
to be completely described by this simple model. In par- 
ticular, the phase with inhomogeneous percolation of FM 
and CO regions is beyond the scope of this paper. 
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